Remarks on Duality Transformations and Generalized Stabilizer States 
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We consider the transformation of Hamilton operators under various sets of quantum operations 
acting simultaneously on all adjacent pairs of particles. We find mappings between Hamilton op- 
erators analogous to duality transformations as well as exact characterizations of ground states 
employing non-Hermitean eigenvalue equations and use this to motivate a generalization of the sta- 
bilizer formalism to non-Hermitean operators. The resulting class of states is larger than that of 
standard stabilizer states and allows for example for continuous variation of local entropies rather 
than the discrete values taken on stabilizer states and the exact description of certain ground states 
of Hamilton operators. 

PACS numbers: 03.67.Hk,03.65.Ud 



I. INTRODUCTION 

The investigation of quantum many-body systems with 
the tools of entanglement theory [Tj has recently re- 
ceived considerable attention addressing long-standing 
questions in the latter employing new methods developed 
in the former. This includes for example the scaling of 
block entropies [|J and geometric entropies a new 
improved understanding and generalization of numerical 
methods such as DMRG on the basis of matrix-product 
states [3, Q and the development of novel approaches 
based on new classes of efficiently describable quantum 
states such as weighted graph states [y, 0] . Most of these 
methods are based on the efficient description of quantum 
states and are taking place in the Schrodinger picture 
while the numerical algorithm described in [6( is a mix- 
ture of Schrodinger and Heisenberg picture, considering 
both transformations of Hamiltonians and description of 
states. In this note consider only transformations on 
Hamiltonians and operators characterizing states. 



II. DUALITY TRANSFORMATIONS 

In the following I will consider the effect of sequences 
of one- and two-qubit quantum gates, such as control 
phase gates and single qubit operations, in a spirit not 
dissimilar to the action of a cellular automaton Q on 
Hamilton operators. Here and in the following we employ 
the notation 



X = 



1 

1 



Y = 







z = 





-I 



(1) 



for the Pauli-operators and I for the identity to simplify 
notation. 

In general, it cannot be expected that with a few sim- 
ple steps of the nature presented in Fig. [T]one can diago- 
nalize a Hamilton operator (even though later on in this 
section such an example will be discussed which is then 



used to motivate the definition of generalized stabilizer 
states). It may however be possible to transform Hamil- 
ton operators into each another, thus establishing useful 
equivalences between seemingly different systems. With 
this aim in mind we would like to explore what effects the 
application of sequences of quantum gates may have on 
Hamilton operators describing quantum many body sys- 
tems. To avoid problems with operator ordering we will 
consider finite translation invariant systems with open 
boundary conditions. This will allow me to use settings 
such as those in Fig. [1] instead of the standard cellular 
automaton in the Margolus partitioning (see e.g. Fig. [5] 
and Q for an excellent exposition of quantum cellular au- 
tomata and their rigorous definition and classification). 
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FIG. 1: Time progresses from bottom to the top. Sequence 
of application of CNOT gates to a set of TV qubits followed by 
the application of Hadamard gates implementing the mapping 
H|0) = (|0> + |l»/\/2 and H\l) = (|0) - |1})/V2. 



The effect of the sequence given in Fig [T] is a transfer- 
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behaviour and strong finite size effects .13] 

, is defined 

L-1 L 

H = J2 JXk-iZ k X k+1 + ^BZ k . (4) 



k=2 



k=l 



where Zj, + i = ILl if the system consists of L spins and we 
use the abbreviation X k = i®^" 1 ) ® X <g> . Con- 

noisseurs will recognize this transformation as a special 
example of a duality transformation Q which can map 
for example the strong coupling regime of a Hamiltonian 
onto the weak coupling regime. Generalizations of such 
duality transformations for Ising models with multi-qubit 
interactions [f Oj may be constructed easily by similar se- 
quences of CNOT gates. Let us now apply this trans- 
formation to the Ising model in transverse magnetic field 
(whose strength we have set to B = 1 for convenience) 
whose Hamiltonian, with open boundary conditions, is 
given by 



Firstly, employing a sequence of controlled phase gates 
between any neighboring pair reveals that the model is 
self-dual in the sense of the Ising model discussed above. 
Employing again the gate sequence in Fig. [T]we find that 
the Hamiltonian eq. (U) is mapped onto 



H 7 



L-1 

E 

k=2 



L-1 



T = BX L - JYkYk+i + BX kXk+i 



(5) 



fe=i 



L-1 L 

H(J) = J2J X kX k +i + Y,Z k . 

k=l fe=l 

Then the above transformation leads to 

L L-1 



(2) 



Again, in the limit L — > oo we expect the correspondence 
between the cluster Hamiltonian and the anisotropic XY 
model to become exact and thus relating the critical be- 
haviour of the two models. Note however that the BXl 
term cannot be neglected in the extreme limit J = 0. In- 
deed, for J = the Hamiltonian in eq. (U|) has a unique 
ground state while eq. (0 without the term BXl term 
would be two-fold degenerate. 

It is also straightforward to see that by the same trans- 
formation the Hamiltonian 



L-1 



L-1 



H T (J) = TH(J)T^ = X L -JZ 1 + Y JZ k + Y X k X k+1 H= Y< 'hX k ^Z k X k+l + ]T J 2 X k X k+1 +Y^BZ k . 

' * ' * k — 2 k—1 = 1 



k=l k=l 



= X L - JZi + JHiJ- 1 ). 



k=2 

is mapped to 



Evidently, the two Hamilton operators are not identical 
but it is reasonable to expect that in the limit L —> oo 
the two terms X^ and JZ\ may be neglected so that we 
can state 



L-2 

E 

k=l 



L-1 



(6) 



BX L - £ JiY k Y k+1 + BX k X k+1 + £ J 2 Z k 



fc=i 



k=2 



H T {J) = JHiJ- 1 ) 



(3) 



Because the operators on both sides obey the same al- 
gebra this may then be viewed as a statement about the 
symmetry of the Hamilton operator itself so that we ex- 
pect for the energy eigenvalues E(J) = Ji?(J _1 ) Q. 
With this relation we may conclude that if the spectral 
gap of the Hamiltonian vanishes for one set of values J 
then it will also vanish for J -1 Thus the assumption 
of a single critical point would then force the conclusion 
that it must be found at J = 1. This turns out to be 
correct in this case and indeed, near the critical point, 
the gap above the ground state for example is given by 
AE = 2 1 1 - J -1 1 thus satisfying eq. ©. The validity of 
these arguments are however perhaps not quite as non- 
trivial as one may expect as it is not per-se clear that a 
small perturbation at the boundaries leaves the spectrum 
essentially unaffected. Such a property requires proof (in 
fact and perhaps not too surprising this is not correct in 
the extreme cases J — and J = oo) even though one 
can expect it to hold in most physically reasonable cases. 

Now we move beyond self-duality and link the cluster 
Hamiltonian [T2, [H[ to the an-isotropic XY model. The 
cluster Hamiltonian, exhibiting some interesting critical 



(7) 

which is the XK-model in a transverse field whose criti- 
cal behaviour is well known. Employing operations such 
as those in figures [T] and [2] one may obtain a large number 
of relationships between Hamilton operators that become 
exact in the asymptotic limit. There are various possible 
directions in which to extend such an approach. One may 
for example consider the transformations that emerge 
from the Trotter decomposition of the time-evolution of 
Hamiltonians H that one has decomposed into two parts 
fi, and H 2 such that H = H x + H 2 and [H X ,H 2 ] = 
[14| . Strictly, speaking it is not necessary to consider 
only unitary operations. Needless to say that in this case 
the spectra of the Hamiltonians are not connected in a 
very transparent way. The next example demonstrates 
that such an approach may nevertheless be useful. The 
following example will also serve to demonstrate that in 
some cases the above approach actually allows in a sim- 
ple way to obtain the exact solution of certain Hamilton 
operators. An example for that, which also serves to mo- 
tivate the definition of generalized stabilizer states, is 
Lemma I — The ground state of the translation in- 
variant Hamiltonian 



N 



H — [— JZ k -\X k Z k+ i + BZ k ] 



(8) 



fc=i 
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with periodic boundary conditions is uniquely determined 
by the N eigenvalue equations of non-Hermitean opera- 
tors 



fe-i 



A 
A" 1 



for all i and 




(9) 



(10) 



Proof: The strategy in the following will be to map the 
Hamiltonian eq. ([5]) to a sum of single particle Hamilto- 
nians that can be solved trivially. To this end, consider 
the operator 



T 



N 

n 



x 1 / 2 

A-V2 



N 



fc+i 



in) 



j fe=i 



where we define the controlled phase Uk,k+i acting on 
the qubits k and k + 1 as 



U, 



k,k+i 



10 
10 
10 
-1 



(12) 



From the observation that the well-known cluster-state 
\Cluster) = JlfcLi Uk,k+i\+)® N is uniquely determined 
as the state satisfying for all i the eigenvalue equations 
Zi-iXiZi + i\Cluster) — [Cluster) [151 ] we immediately 
find that 



I*) = T\ 



(13) 



satisfies the eigenvalue equations ©. The transformed 
Hamilton operator Ht = T~ 1 HT is given by 



N 

fc=i 



B 

-JX 



-JA- 1 
—B 



(14) 



Now we need to determine A such that the ground state 



of Ht is given by n,=i This is easily found to be 



and 



(15) 




E = -Ny/B 2 + J 2 . 

The observation 

H\V) = E \V) & H T \+)® N = E \+)® N . 
then completes the proof. 



(16) 
(17) 



Remark I - We could have also applied a different 
transformation R which applies the local and non-local 
operations in reverse order to that in T. Indeed 



N 



N 



R=n Uk.k+i n Uj 

3=1 



(18) 



fc=i 



where Uj is the unitary transformation that makes 



U] 



B -J 

-J -B 



U, 



(19) 



diagonal and smaller eigenvalue being the vector |1) = 
( , ) . The transformation on the Hamilton operator is 
presented pictorially in Fig. O Thus the state \$f) of 
eq. © may also be written as |^) = R\l). It is note- 
worthy that by interchanging the order in which the sin- 
gle particle and the two-particle operations are applied 
changes the single particle operator from a non-unitary 
to unitary. This approach also allows immediately for 
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FIG. 2: Time is progressing from bottom to the top. A cel- 
lular automaton first applies controlled phase gates (they are 
symmetric with respect to interchange of control and target) 
between nearest neighbours and then single qubit gates U 
such that the combined action maps the Hamilton operator 
eq. |H} onto a single particle Hamiltonian that is diagonal in 
the computational basis. 

the determination of the entire eigenvalue spectrum. The 
energy separation between eigenstates is simply given by 
the splitting between the eigenvalues of [B — J; — J — B]. 
Thus it becomes evident that the Hamilton operator does 
not possess a quantum critical point. 
Remark II - For J = the ground state is a product 
state, while for B = the ground state is a cluster state. 
The transition between those two regimes is continuous 
and thus we observe that the states |^) can realize local 
entropies ranging continuously in the interval [0, 1]. This 
is in contrast to standard stabilizer states whose local 
entropies are quantized in integer units which in itself 
is enough to see that standard stabilizer states will be 
ground states to a very limited set of Hamiltonians (see 
|24| for a much more detailed discussion). 
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The simple example above suggests that the exten- 
sion of the stabilizer formalism to non-Hermitcan op- 
erators may be useful. To remind the reader, a 'stabi- 
lizer operator' for an N-particle state is a tensor product 
of N operators picked from the set {X, Y, Z,l}. A set 
G = {gi, . . . , <7at} of ./V mutually commuting and inde- 
pendent stabilizer operators 16] is then called a 'gener- 
ator set' uniquely identifying a state |^) that satisfies 
<7fc |\&) = |W) for all k — 1,...,N. Now we generalize 
these notions to the definition of generalized stabilizer 
states. 

Definition I (Generalized stabilizer states) - For 

N qubits, a generalized stabilizer state |W) is the unique 
eigenstate to eigenvalue 1 of the N mutually commut- 
ing and independent 'generalized stabilizer operators ' 
{<?i, . . . , gw} where each g^ is an N-fold tensor product 
of arbitrary, possibly non-Hermitean, linear operators. 

From this definition we can immediately draw some 
straightforward conclusions on some general classes of 
states that admit such a generalized stabilizer state de- 
scription. 

Remark III - Note that any pure two-qubit state can be 
written as a generalized stabilizer state. This follows di- 
rectly from the fact that each pure two qubit state can 
be obtained from a maximally entangled singlet state by 
the local application of linear operators. For three qubits 
all members of the GHZ class are generalized stabilizer 
states while the members of the W-class may be approx- 
imated arbitrarily well by generalized stabilizer states 
(this follows from the well-known classification of pure 
three-qubit states under the action of local linear maps 
Some interesting results have been obtained re- 
cently in [IliHiiHH concerning the entanglement content 
of certain stabilizer states and it would be interesting to 
see how these results may be generalized to the setting 
of generalized stabilizer states. 

Note that one may extend the class of states which 
admit a description also in another direction, namely by 
grouping together neighbouring qubits. For kN qubits, a 
generalized stabilizer state |^) is the unique eigenstate to 
eigenvalue 1 of the N mutually commuting and indepen- 
dent 'generalized stabilizer operators' {g\, . . . , giy} where 
each is an N-fold tensor product of arbitrary, possibly 
non-Hermitean, linear operators acting on k qubits. It 
is evident that this will, for k = 3 for example, describe 
all possible three-qubit states but of course this comes at 
the expense of an exponential increase in complexity of 
description. In some sense the above concepts of gener- 
alized stabilizer states are included in other descriptions 
of quantum states. In fact, for qubits deformed weighted 
graph states as they have been introduced in [|| (see also 
uL [22| for the concept of weighted graph states) incorpo- 
rate the generalized stabilizer states according to Defini- 
tion I. However, in this picture the description is again on 
the level wave-functions (23j |, i.e. the Schrodinger picture, 
while the present approach, via eigenvalue equations, is 



situated in the Heisenberg picture. The weighted graph 
state picture of @, 0, [HJ is probably too general to al- 
low for a detailed quantification of multi-particle entan- 
glement while the generalized stabilizer states may well 
admit more detailed results for admittedly a less gen- 
eral class of states. This will be the subject of a future 
publication. 

As a final remark it should be noted that by definition, 
generalized stabilizer states possess a unique characteri- 
zation of the quantum state of an n-qubit system employ- 
ing only resources that are polynomial in n. This alone, 
however, is not sufficient for applications. It is also im- 
portant to be able to derive relevant physical quantities 
directly from the stabilizer formalism. Indeed, having 
first to deduce the state explicitly and then computing 
the property from the state would generally involve an 
undesirable exponential overhead in resources. While one 
can expect a direct approach to be possible in principle, 
it is evident that detailed and explicit presentations of al- 
gorithms to achieve these tasks in a systematic way and 
whose convergence is proven are of interest. For standard 
stabilizer states such a program has been carried out and 
detailed algorithms have been provided [HI]. There, nor- 
mal forms in the context of bi-partite entanglement have 
been found together with algorithms with proven con- 
vergence to obtain these. These tools may be transferred 
directly to generalized stabilizer states whenever they are 
obtained from standard stabilizer states by the local ap- 
plication of linear operators. 



III. CONCLUSIONS 

In this work we have studied duality relations between 
Hamilton operators from the viewpoint of sequences of 
quantum operations. The close relationship of these 
transformations to duality transformations in Hamilton 
operators has been noted. The same approach has then 
been used to derive an exact solution for a many-body 
Hamiltonian which has led to the characterization 
of the ground state via eigenvalue equations using 
non-Hermitean operators. This motivated the definition 
of a generalization of the concept of stabilizer states to 
incorporate non-Hermitean operators. The properties of 
such states and the potential offered by this approach 
still remain to be explored. 
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